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Abstract 



We compute the ratio A^/Ajjg, where the scale parameter A^ is associ- 
ated with a lattice formulation of QCD. We consider a 3-parameter family 
of gluon actions, which are most frequently used for 0(a) improvement a la 
Symanzik. The gluon action is put togeter with standard discretizations for 
fermions (Wilson/clover, overlap), to provide for several possible combina- 
tions of fermion and gluon actions. We employ the background field technique 
in order to calculate the 1PI 2-point function of the background field; this 
leads to the coupling constant renormalization function, Z g , at 1-loop level. 

Our results are obtained for an extensive range of values for the Symanzik 
coefficients. 
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I. INTRODUCTION 



The A parameter of QCD has been a subject of interest for almost three decades, since it is 
the necessary "yardstick" needed to convert dimensionless quantities coming from numerical 
simulations into measurable predictions for physical observables. 

Ever since improved gluon and fermion actions started being employed more frequently 
in numerical simulations, a number of calculations of the A parameter on the lattice have 
been carried out, using various techniques and discretization prescriptions. Older results 
involving Wilson gluons [1], Wilson/clover fermions [2,3], overlap fermions [4] can be found 
in the literature. Some recent results regarding domain wall fermions can be found in Ref. 
[5]- 

A calculation of the A parameter which is missing is the one involving the Symanzik 
improved gluon actions which are widely used in recent simulations. The task of the present 
work is to fill this gap, while at the same time we confirm some of the existing results 
mentioned before. The contribution of fermions to the computation at hand is independent 
of the choice of gluon action; similarly, gluon contributions do not depend on the fermion 
action. This fact will enable us to combine our results with previous findings regarding 
Wilson/clover fermions [3] and overlap fermions [4]. 

The scale parameter, A L , associated with a lattice formulation of QCD provides a relation 
between the lattice spacing, a, and the bare coupling constant g Q . It is a particular solution 
of the renormalization group equation, taking the form 



qAl = exp 



-/ 



9° dg 



exp 



2\-b 1 /2b 



1 + 0{gl) 



(1) 



where Pi(go) is the lattice /9-function, and bo, b\ the first two coefficients of its perturbative 
expansion 



bo 



11 2 

<*0» KT N ' ~ 3 N > 
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(2) 
(3) 



(Nf-. the number of fermion flavors, N c : the number of colors.) 

The A parameter is a dimensionful quantity; as such it cannot be directly obtained from 
the lattice. Instead, the quantity which is calculable is the ratio between A L and the scale 
parameter in some continuum renormalization scheme such as MS: Al/Aj^. To this end, 
it suffices to compute the coupling constant renormalization function Z g , relating the bare 
lattice coupling g Q to the MS'-renormalized coupling g. 

For the purposes of this calculation, we employ the background field technique [6-8]. 
This technique lends itself particularly well to evaluating Z g , since it obviates the need to 
consider any 3-point functions. 

The rest of the paper is organized as follows: In section II we present all the necessary 
background and set up. Our results are shown in section III and, finally, a brief discussion 
regarding some aspects of our calculation and findings is contained in section IV. 
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II. FORMULATION OF THE PROBLEM 



We use the Symanzik improved gauge field action, involving Wilson loops with 4 and 6 
links 1 . In standard notation, it reads [9] 



g 2 



c £ Re Tr (1 - C/ plaq ) + c x £ Re Tr (1 - C/ rect ) 

plaq 



rect 



+ c 2 £ Re Tr (1 - C/ Chair ) + c 3 £ Re Tr (1 - C/ paral ) 

chair paral 



(4) 



The lowest order expansion of this action (together with the gauge fixing term, with 
gauge parameter £, see Eq.(lO)), leading to the gluon propagator, is [10] 



^ 0) = 



i ri* d A k 



w/a (2?r) z 



/if 



(5) 



where: 
and: 



G fiu (k) — kpk v + {k 2 p 5^ v kpjz p 8 pi ^ d^p 



d 



fW v l - V) [C - d a 2 P - C 2 a 2 (A; 2 + A; 2 ) 

- > = : 



k 



■ sm 



a 



k 2 — 



The coefficients Cj are related to the Symanzik coefficients q by 

C = c + 8ci + 16c 2 + 8c 3 , Ci = c 2 + c 3 , C 2 = c\ - c 2 - c 3 



(6) 



The Symanzik coefficients must satisfy: C = 1, in order to reach the correct classical 
continuum limit. 

Regarding the fermion part of the action, a variety of discretizations are presently used 
in Monte Carlo simulations. The contribution of fermions to 1 loop is independent of the 
regularization chosen for the gluonic part; vice versa, gluon contributions do not depend 
on the fermion action. Consequently, the results of the present work can be directly com- 
bined with those of previous calculations regarding Wilson/clover fermions [3] and overlap 
fermions [4], yielding the A ratio for a variety of possible combinations of fermion and gluon 
actions. 

In the background field method, link variables are decomposed as [6] 

U„{x) = V^x) U c ,(x) (7) 
in terms of links for a quantum field and a classical background field, respectively 



1 1 x 1 plaquettes, 1x2 rectangles, 1x2 chairs (bent rectangles), and parallelograms wrapped 
around elementary cubes. 
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V^x) = e i9Q ^ x \ U Cf ,(x) = e iaB ^ x) (8) 
The N c x N c Hermitian matrices and can be expressed as 

Q,(x)=t a Q;(x), B^x)=fB*(x), Tr[t a t b ]= l -5 ab (9) 



A choice of gauge is required for the perturbative expansion; an appropriate gauge-fixing 
term is 

Sgf = TZ7 E E Tr^Q^Q,] (10) 

This term breaks gauge invariance with respect to Q^, as it should, but succeeds in keep- 
ing the path integral as a gauge invariant functional of B^. The definition of the lattice 
derivative, which is covariant with respect to background gauge transformations, is 

D~(U c )Q u {x) = U~^(x - ey)Q v {x - e^)U c ^(x - e M ) - Q v (x) (11) 

Since the quantities we will be studying are gauge independent, we chose, for convenience, 
to work in the Feynman gauge, £ = 0. Covariant gauge fixing produces the following action 
for the ghost field oo 

s 9h = 2Y J Y.^( D+ M x )) ] ( D+ M x ) + % 9* [QMMx)\ + lm [q„(x),d+u(x) 

X fJL 

-^9 2 [Q,(x),[Q,(x),D;u;(x)]]+---), (12) 

where D+u(x) = U CjJi {x)uj{x + fi)U~^{x) — uj{x). 

Finally the change of integration variables from links to vector fields yields a Jacobian 
that can be rewritten as the usual measure term S m in the action 

^ = E{^Tr(Q M (x) 2 ) + ") (13) 

The measure part will not contribute to the present calculation. 

In order to compute we need to evaluate the renormalization function Z g for the 
coupling constant, up to 1 loop 

9o = Z g (g ,afi)g (14) 



where ji is the renormalization scale in the MS scheme. Writing 

Z g ( 9o , ajif = 1 + g 2 o (2b \n(ajl) + l Q ) + 0(g 4 ) (15) 

one has 

/ o = 26 ln(A L /A 1 ^) (16) 
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To obtain Z g we only need to calculate the one-particle irreducible (1PI) 2-point function 
of the background field, r( 2,0 '°)(p, — p)^ u , on the lattice, to one loop. Color symmetry and 
lattice rotational invariance allow one to write [11] 



E r (2A0) (P, -v)t = - 3 ^ 2 [i - "(p)1 1st 



(17) 



where u{p) is a Lorentz invariant amplitude on the lattice, up to terms which vanish as 
a — > 0; v{p) is perturbatively expanded as 



^(p) = E^ w (p) 



(18) 



i=l 



The background field formalism has the advantage that Z s is directly related to the 
background field renormalization function Z A , through: Z 9 (g Q , a/l) 2 Z^(g , a/i) = 1- Conse- 
quently, no 3-point functions are needed for the evaluation of Z g . In terms of v(p), one can 
express Z g as 



Z g ( 9o , a^f = l + gl (4\p/fi) ~ + O(gt) 



where 



16tt 2 



1 



11 p 2 205 3 
T n ^ + ^6 r + 2(r^O + 4(l-0 2 



+ 



16tt 2 



2 

3 ln ^ 
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(19) 



(20) 



is the analogous 1-loop amplitude in the MS scheme. 



III. COMPUTATION AND RESULTS 

The Feynman diagrams shown in Fig. 1 contribute to v^\p). All algebraic manipulation 
of these diagrams was performed automatically using our software written in Mathematica. 
Once we have computed u^(p), we use Eqs.(15) and (19) in order to obtain l Q , which is the 
sum of a part involving only the gluon and ghost action, and a part involving the fermion 
action, i.e., 

l = 11 + N, ■ ll (21) 




Fig. 1. One-loop diagrams contributing to r^ 2 ' 0,0 ). A wavy (solid, dashed) line 
represents gluons (fermions, ghosts). The letter B stands for the external background field. 
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The dependence of 1% on the Symanzik coefficients is rather complicated and cannot 
be given in closed form. However, given that the gluon propagator depends only on the 
combinations Ci, C 2 (c.f. Eqs.(4), (6)) we can reexpress all diagrams in terms of Ci, C 2 
and one additional parameter, say, c 2 ; in this case the dependence on c 2 (at fixed Ci, C 2 ) is 
polynomial. Thus, the part of / Q involving gluon and ghost fields, 1%, can be written as 

1 c 

l 9 = ai — + a 2 N c + a 3 + a 4 c 2 N c + a 5 c^ c (22) 

where aj are numerical constants (dependent on Ci, C 2 ) evaluated via numerical integration 
over loop momenta. We consider ten sets of different values for the Symanzik coefficients, 
corresponding to the most commonly used actions, shown in Table I: The plaquette ac- 
tion, the tree-level Symanzik improved action, the Luscher-Weisz tadpole improved actions 
(T1LW), the Iwasaki action and the DBW2 action (see [12-17]). The quantities Oj, for each 
one of the ten sets of parameters, are presented in Table II. The variable c 2 can be freely 
varied; c , C\ and c 3 are then adjusted accordingly so as to keep C , C\ and C 2 fixed. 

The fermionic part of / Q , denoted by l[, was calculated in [4] using Neuberger's overlap 
formulation of chiral fermions [18], leading to 

« = " MP) (23) 

where kf(p) varies from 0.07 to 0.08 in a typical range of the overlap parameter p. For an 
extended list of values of kf(p) see Table I of Ref. [4]. 

In order to assess quantitatively the effect of the Symanzik improved actions on the A 
parameter, one may consider the ratio 



/ a / a \Symanzik r 

l A L/ A M57 _ J_ fjSymanzik _ ,Wilson\ 



(24) 



This quantity is independent of the fermion action but still depends on the number of flavors, 
Nf, through b . For completeness, we report the value of l Q found in the literature for Wilson 
gluons and Wilson/clover fermions [1,3] 

l = 1/(SN C ) - 0.169955999 N c + N f l 01 (25) 



where : l m = 0.006696001(5) - 0.00504671402(1) c sw + 0.02984346720(1) c 2 sw (26) 

In the above, the Wilson parameter is set to r = 1 and the clover parameter, csw, can be 
chosen arbitrarily; the dependence on csw is seen to be polynomial. 

In Table III we list the values of the ratio r A for Nf = and Nf = 2. We present r A 
for each set of parameters shown in Table I, setting N c = 3 and c 2 = 0. We also list the 
A ratio, A L / Ajjg\ Nf=0 , in the pure gauge theory, and with two flavors of Wilson fermions, 
Ai/Ajg^-|7v /=2 . We stress that r^(Nf = 2) is the same for all types of fermion actions. 

In Fig. 2 we present our results for 1% as a function of both C\ and C 2 , for N c = 3 and 
c 2 = 0. The range of values for C\ and C 2 was selected so as to encompass all values used in 
current simulations. We can see that the dependence on C\ is almost linear while dependence 
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on C 2 is more complicated. The crosses correspond to the ten actions shown in Table I. In 
Fig. 3 we plot the ratio defined in Eq.(24) as a function of C 2 . Once again, we have set 
N c = 3, c 2 = c 3 = and thus C 2 = C\. 




Fig. 2. l as a function of the parameters C\ and C 2 (N c = 3, c 2 = 0). The crosses 
denote the 10 set of parameters, identified by their C\, C2 values, as shown in Table I. 
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Fig. 3. ta as a function of C 2 , for Nj = 0, 2. 
We have set c 2 = C3 = (and therefore C 2 = ci) and N c 
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For easier comparison we report some existing numbers for the ratio Al/Ajj^ using 
Wilson gluons and Wilson or overlap fermions (see e.g. Refs. [3,4]) 

Wilson gluons, Wilson fermions : Al/Ajj^ = 0.0243589 

Wilson gluons, overlap fermions (p = 1.0) : Al/Aj^ = 0.0172702 

Wilson gluons, overlap fermions (p = 1.4) : Al/Aj^ = 0.0172317 



IV. DISCUSSION 

In the present work we evaluated Al/Aj^, for a 3-parameter family of Symanzik im- 
proved gluon actions; to this end, we computed Z g , up to 1 loop, using the background 
field technique. Only diagrams with two external background fields, corresponding to the 
1PI two-point function of the background field, were involved in the calculation, as shown in 
Fig. 1. Alternatively, one could study Z g by considering the gluon-gluon-gluon, gluon-ghost- 
antighost or gluon-fermion-antifermion three-point functions, together with the self-energy 
diagrams for the gluon, ghost and fermion fields; of course, the computation would be much 
more cumbersome in this case, due to the complexity of the Symanzik improved actions, 
resulting in lengthier algebraic expressions. It is this very fact pointing out the advantage 
of the background field technique. 

All calculations have been performed in the Feynman gauge (£ = 0), and the conversion 
of lengthy integrands (~100,000 terms) into an efficient Fortran code was carried out by our 
"integrator" program, a metacode written in Mathematica. The numerical integrals were 
evaluated for lattices up to L = 128; the results were then extrapolated to L — > oo. Given 
that only a restricted set of functional forms is sufficient to describe the behavior of the 
results with L, the systematic error resulting from such an extrapolation can be estimated 
quite accurately. 

Special attention was given to the extraction of the dependence on the external momen- 
tum p. The algebraic expressions coming from the evaluation of Feynman diagrams were 
split into a logarithmically divergent part, comprised of a limited set of tabulated lattice 
integrals, and a (much larger) part which is Taylor expandable up to second order in p. We 
have seen explicitly that terms of order O(p ) cancel upon summation of gluon, ghost and 
fermion diagrams separately, compatibly with gauge invariance. 

Our results are functions of the Symanzik coefficients C±, C 2 and c 2 . At fixed C\ and C 2 , 
the dependence on c 2 is seen to be a second order polynomial, thus no particular values of c 2 
have to be chosen a priori; conversely, to investigate the effect of the remaining coefficients, 
we selected a mesh of 25x27 values of C±, C 2 for numerical integration. The dependence on 
C\ turns out to be almost linear, while the C 2 dependence is more complicated (see Fig. 2). 

Given that the gluon and fermion parts of the action give disjoint contributions to l D , 
our present result can be directly combined with contributions from a variety of different 
fermion actions, to yield the complete effect on A L . The number of colors, N c , and the 
number of fermion flavors, Nf can be chosen arbitrarily. 

Through Eq.(24) one can assess the effect of the Symanzik improved actions on the A 
parameter. All of the actions shown in Table. I, with the exception of the DBW2 action, give 
similar results (of order 10 1 ) for r\. The most drastic effect on the A parameter originates 
from the DBW2 improved action, where r A is of order 10 3 . 
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TABLES 



TABLE I. The coefficients cq, c\, C3 (02 = 0), corresponding to some of the most 
commonly used actions, along with the respective values for C±, Ci- 







Action 




Co 


Cl 




C3 




Ci 


c 2 


Set 


1: 


Plaquette 




1.0 


0.0 


0. 


.0 


0. 


.0 


0.0 


Set 


2: 


Symanzik 




1.6666666 


-0.083333 


0. 


.0 


0. 


.0 


-0.083333 


Set 


3: 


TILW, (3 = 


8.60 


2.3168064 


-0.151791 


-0. 


.0128098 


-0. 


0128098 


-0.138981 


Set 


4: 


TILW, = 


8.45 


2.3460240 


-0.154846 


-0. 


0134070 


-0. 


0134070 


-0.141439 


Set 


5: 


TILW, = 


8.30 


2.3869776 


-0.159128 


-0. 


.0142442 


-0. 


.0142442 


-0.144884 


Set 


6: 


TILW, p = 


8.20 


2.4127840 


-0.161827 


-0. 


.0147710 


-0. 


.0147710 


-0.147056 


Set 


7: 


TILW, p = 


8.10 


2.4465400 


-0.165353 


-0. 


.0154645 


-0. 


0154645 


-0.149889 


Set 


8: 


TILW, p = 


8.00 


2.4891712 


-0.169805 


-0. 


.0163414 


-0. 


.0163414 


-0.153464 


Set 


9: 


Iwasaki 




3.648 


-0.331 


0. 


.0 


0. 


.0 


-0.331 


Set 


10: DBW2 




12.2688 


-1.4086 


0. 


.0 


0. 


.0 


-1.4086 



TABLE II. Values of the coefficients a\, 02, (13, 04, 05. Set 1 through Set 10 correspond to C\, C2 
shown in Table I. 



Set 



ai 



0-2 



03 



04 



(Ik 



1 0.12499999997(6) 

2 0.04217165191(7) 

3 -0.0082581838(1) 

4 -0.010122622374(4) 

5 -0.01268965657(2) 

6 -0.0142802781(2) 

7 -0.01632979450(5) 

8 -0.01886971906(3) 

9 -0.07528696825(4) 

10 -0.204424737(1) 



-0.1699559990(1) 
-0.0833756545(3) 
-0.0310360175(6) 
-0.0290818603(4) 
-0.0263884374(3) 
-0.0247177255(2) 
-0.0225635251(5) 
-0.01989088288(9) 
0.0433593330(5) 
0.19876966(7) 



0.43112525414(6) 

0.31095652446(7) 

0.24637322205(3) 

0.244114928978(2) 

0.241020561259(1) 

0.23911193181(8) 

0.2366619590(1) 

0.23364084600(6) 

0.17401920011(2) 

0.06102725834(6) 



-0.0958290656(4) 
-0.031584124(1) 
0.000854777(6) 
0.000128847(6) 
0.001465153(3) 
0.0022825409(6) 
0.003324978(1) 
0.004598716(5) 
0.021865609(1) 
0.03791059(2) 



-0.6576721162(8) 

-0.402576126(2) 

-0.283219286(3) 

-0.279256305(1) 

-0.2738512399(6) 

-0.270531792(3) 

-0.2662885752(6) 

-0.2610823949(1) 

-0.159911864(1) 

-0.027913991(2) 
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TABLE III. The ratio r\ defined in Eq.(24) for each set of parameters and for Nf = 0, 2, 
along with the respective values for A^/ Aj^g\j^ f= Q and Ai J /Ajrjs\N f =2- We have set N c = 3 and 
c 2 = 0. 







A rt i on 




r\ (Nt — 0) 


A r / A 1 at n 


1 A y 1 v / — 


A r / A 1 AT o 


Set 


1: 


Flaquette 




1.00000 


0.034711 


1.00000 


0.024359 


Set 


2: 


Symanzik 




5.29210 


0.18369 


6.65946 


0.16222 


Set 


3: 


TILW, p = 


8.60 


14.4779 


0.50254 


20.9316 


0.50987 


Set 


4: 


TILW, = 


8.45 


15.0329 


0.52181 


21.8471 


0.53217 


Set 


5: 


TILW, = 


8.30 


15.8330 


0.54958 


23.1751 


0.56452 


Set 


6: 


TILW, P = 


8.20 


16.3507 


0.56755 


24.0392 


0.58557 


Set 


7: 


TILW, P = 


8.10 


17.0432 


0.59159 


25.2012 


0.61387 


Set 


8: 


TILW, p = 


8.00 


17.9435 


0.62284 


26.7214 


0.65090 


Set 


9: 


Iwasaki 




61.2064 


2.1245 


107.957 


2.6297 


Set 


10: DBW2 




1276.44 


44.306 


3423.05 


83.382 
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